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PREFACE 


The  research  discussed  in  this  report  was  sponsored  by  the  Defense 
Advanced  Research  Projects  Agency  and  is  intended  to  provide  a simple 
analytical  means  of  estimating  the  effects  of  wall  boundary  conditions 
on  boundary-layer  stability.  Accuracy  was  not  the  primary  concern  of 
the  study.  However,  emphasis  was  placed  on  isolatin;^.  the  critical 
boundary-layer  parameters  and  providing  some  quantitative  prediction 
of  their  effects  on  boundary-layer  stability. 
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SUMMARY 


An  integral  solution  of  the  two-dimensional  boundary-layer  equa- 
tions for  water  with  pressure  gradient,  heat  transfer,  and  suction  was 
used  to  investigate  laminar  boundary-layer  stability.  There  are  theo- 
retical bases  to  indicate  that  the  effects  of  suction,  wall  heating, 
and  p~essure  gradient  on  critical  Reynolds  numbers  could  be  correlated 
as  a function  of  universal  bound ary- layer  parameters.  Either  the 
boundary- layer  shape  factor,  H,  or  the  parameter  proposed  by  Wazzan, 
u"(0),  could  be  used.  These  two  parameters  are  not  Independent,  The 
preference  of  one  parameter  over  the  other  depends  on  the  nature  of 
the  data  to  be  correlated.  There  is  apparently  no  preference  for  theo- 
retical data,  and  H is  preferable  for  experimental  data  because  the 
displacement  and  momentum  thicknesses  are  easier  to  deduce  from  measure 
ments  than  the  wall-shear  parameters. 

The  results  of  the  analysis  are  used  for  a qualitative  discussion 
of  the  effects  of  suction,  wall  heating,  and  pressure  gradient  on 
boundary- layer  stability.  Suction  is  the  most  effective,  and  pressure 
gradient  the  least  effective,  means  to  stabilize  a boundary  layer. 
Further,  the  stabilizing  effects  of  suction  and  wall  heating  are  en- 
hanced by  a favorable  pressure  giadient.  The  application  of  this 
analysis  to  calculations  of  the  critical  Reynolds  number  for  an  arbi- 
trary body  shape  is  disc_ssed. 


f 
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SYMBOLS 


a = velocity  profile  parameter,  Eq.  (10) 
b = temperature  profile  parameter,  Eq.  (11) 
c = wave  speed 
c = specific  heat 

= boundary-layer  shape  factor,  Eq.  (7) 

P = pressure 
Pr  = Prandtl  number 


Q = suction  parameter,  Eq.  (12) 
I = rrii-ira1  Revnolds  number 


c 

T = temperature 

u = tangential  velocity  in  the  boundary  layer 

U = tangential  velocity  external  to  the  boundary  layer 

V = normal  \elocity  in  the  boundary  layer 

V = wall  suction  velocity 
w 

X = tangential  coordinate 
y = normal  coordinate 
a = wave  number 
6 = momentum  layer  thickness 
6 = displacement  thickness,  Eq.  (7) 

6^  = thermal  layer  thickness 
A = 6t/6 

E = variable  viscosity  parameter,  Eq.  (12) 
p = molecular  viscosity 
A = Pohlhausen,  shape  factor,  Eq,  (12; 

K = thermal  conductivity 
p = density 

n = boundary  layer  coordinate 
0 = momentum  thickness,  Eq.  (7) 

0 = energy  thickness,  Eq.  (7) 

= variable  viscosity  parameter,  Eq.  (12) 
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Subscrlpts 

e = external  to  the  boundary  layer 
t * thermal  layer 
w * wall 

X = derivative  with  respect  to  tangential  coordinate 
y = derivative  with  respect  to  normal  coordinate 


I . INTRODUCTION 


In  the  study  of  houndary-layer  control,  it  is  well  known  that  suc- 
tion can  completely  stabilize  a laminar  boundary. It  is  also  known 

(2  3) 

that  heating  is  stabilizing  for  a laminar  water  boundary  layer.  ’ 

Both  of  these  effects  have  been  studied  numerically.  But  as  with  most 
numerical  studies,  the  interplay  among  critical  parameters  is  not  illum- 
inated. The  expense  of  conducting  either  a numerical  or  a physical 
experiment  deters  wide  coverage  of  the  various  parameters  that  influ^'ace 
boundary-layer  transition. 

It  is  the  purpose  of  this  study  to  present  a simple  integral  solu- 
tion of  self-similar  flows  that  displays  the  relationship  among  the 
various  boundary-layer  profile  parameters.  These  results  will  sliow 
that  the  critical  Reynolds  number  as  a function  of  pressure  gradient, 
wall  temperature,  and  wall  suction  can  be  correlated  with  a single 
boundary-layer  parameter.  The  assumption  of  local  similarity  will  be 
used  in  a discussion  of  the  calculation  of  the  critical  Reynolds  number 
for  an  arbitrary  body  shape. 
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II.  ANALYSIS 


BOUNDARY-LAYER  FLOW 

In  this  section,  I shall  develop  the  equations  that  define  t ae 
boundary-layer  profile  parameters  using  an  integral  analysis.  The 
boundary-layer  stability  parameters  will  be  determined  from  the  approx- 
imate Schlichting-Ulrich  correlation  and  the  Lin  equations  and  the 

profiles  derived  in  this  section. 

Consider  the  flow  of  water  around  a heated  two-dimensional  body. 
The  two-dimensional  laminar  boundary-layer  equations  for  an  incompres- 
sible, variable-viscosity  flow  are: 


5u  ^ . 
3x  ^ Dy 


0 . 


(1) 


where  u and  v are  the  tangential  and  normal  velocities,  x and  y are 
the  coordinates  tangential  and  normal  to  the  body,  the  subscripts  x 
and  y denote  the  derivatives  with  respect  to  x and  y,  "nd  y is  the  var 
iable  viscosity  that  depends  on  the  temperature.  The  specific  heat, 
c , the  density,  p,  and  the  thermal  conductivity,  <,  are  assumed  con- 
stant for  heated  water.  The  appropriate  boundary  conditions  are: 


y = 6:  '1  = Ug,  Uy  - “yy  = T = T^,  Ty  T.^.y  » 0 ; 


^ du  3T  3u  , 3^  . 

0:  u = 0,  T = T^,  pv^Uy  = pU^  ’ 


3T  3^T 

PV  ^ = K ^ . 

" 3y^ 


(4) 


where  v denotes  th»*  wall  suction  velocity, 
w 

The  integral  equations  are  obtained  by  integrating  Eqs.  (2'  and 
(3)  across  the  boundary  layer.  The  results  are: 


M 

d6  'ex^ 
dx  U 

e 


6 (2  + H) 


^w  3u 

„2  9y 

pU  •'w 
e 


(5) 


dx 


T V + 
e w 


K 


(6) 


where 


1 assume  that  the  velocity  and  temperature  profiles  can  be  represented 
by  fourth-order  polynomials  as  expressed  by  the  following: 


^ = a(n  - 3n^  + 3n^  - n^)  + 6n^  - + 


3n 


(8) 
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1 - T 


2 . o_3  „4.  ,2  Q 3 , „ 4 


-~  = b(n^  - 3n;  + 3n^  - n^)  + 6n^  - 8n^  + 3n^  , 


(5) 


e w 


whe 


re  n = y/<5  and  r,^  = y/6^.  The  profil.'!  parameters  are  defined  by 


where 


= X + 1.2 
^ 6 - (I)  + Q ’ 


(10) 


b = 


^ ^ A Pr  ’ 

6 + Q — 


(11) 


pv 


Q = 


/du\  ~ '^w\  6b 

° \dT/  \ y I 6_ 

' /w  \ w / T 


Eb 

A ’ 


(12) 


dU  ^ 6_ 

X = ^ 5^  A = / . 

y dx  o 

w 


For  wall  heating  and  a heated  water  boundary  layer  <})  > 0,  and  for  suc- 
tion Q < 0.  Therefore,  an  increase  in  either  or  both  of  these  param- 
eters produces  an  increase  in  a. 

So  far  the  analytical  procedure  is  similar  tc  that  used  by  Haupt- 
mann, except  that  here  emphasis  is  also  being  placed  on  the  effects 
of  pressure  gradient  and  suction  on  boundary- layer  stability. 

If  the  profiles  given  by  Eqs.  (8)  and  (9)  are  used  to  evaluate  the 
parameters  in  the  boundary-layer  equation,  one  finds  that 


2 

0 4. a_ a 

5 * '35  (5)  (21)  ■ (36)  (7)  ’ 


(13) 


5 * _ 3 

6 


20 


(14) 
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0 a_  (a 

6^  " 1,35  20  ■ (35) 


b)17  ab 

(24)  ■ (35)  (7) 


1 U (T  - T 

J »•  e w 


) Pr  - 1 , (15a) 


« A fpr  (6  - b)  U (T  - T ) 
0 60  e e w 

T 


PR  » 1 


A « 1 


(15b) 


Equations  (15a)  and  (15b)  represent  the  appropriate  cases  for  a heated 
water  boundary  layer.  Substituting  Eqs.  (8)  and  (9)  Into  Eqs.  (5) 
and  (6) 


U 

ex 

U 

e 


0(2  + H) 


(a  + Q)  , 


(16) 


d0 

dx 


(Te  - 


T ) 

w 


K 

pc. 


[b  + AQPr] 


(17) 


For  the  self-similar  Falkaer-Skan  flows,  a,  b,  Q,  and  (J)  are 
constants,  and  x”.  The  wall  temperature  Is  a constant,  and  the 
slight  variations  In  will  be  neglected. 

Solving  Eqs.  (16)  and  (17)  for  6 and  A yields 


2xy  (a  + Q) 
w 


pU^(i  + 4m  t 2mH) 


35  (5)  (21)  (36)  (7) 


(18) 


and 


A 


2 


6U 

e 


2xy^{b  + AQPr} 

2rr4  a (a  -b  b)17  ab  1 
^^[  35  20  (35)  (24)  (36)  (7)J 


Pr  s:  1 (19a) 


3 


pU  6 
e 


Pr[b  + AQPr] 


Pr  » 1 . 


(19b) 
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2 

These  equations  are  not  explicit  expressions  for  6 and  A because  the 
boundary-layer  parameters  a,  b,  and  Q contain  either  i5  or  6^  or  both. 
Of  course,  they  could  be  solved  numerically.  However,  to  find  the 
important  parameters  1 shall  use  simplifying  assumptions  of  local 
similarity. 

With  flat-plate  local  similarity. 


6 


A 


A 


2 

0 


0 


0 


315 

pU^  37 


Pr 


-1/2 


Pr 


-1/3 


Pr  1 


PR  » 1 . 


(20) 


Using  the  above  values  in  the  equations  for  a,  b,  Q,  and  (J),  one 
can  derive  first  approximations  for  these  parameters.  For  example, 
the  case  where  Pr  1 results  in 


^0  = 


e 2Pr 


210m 

•^7 


1/2 


+ 2 


1 - 


4(315) (pU^x) 


37u 


w 


(21) 


O' 

o 

II 

V /4(315)(pUx) 

W f c 

1 + 

U W 37y  I/- 

el  ^w  1/2 

6 

(22) 


The  above  equations  are  the  first  approximation  of  the  parameters 
a and  b.  The  procedure  for  calculating  6,  A,  a,  b,  X,  and  H for  arbi- 
trary body  shape  and  arbitrarily  distributed  but  modest  mt^gnitude  of 
suction  is  to  calculate  a and  b from  Eqs.  (21)  and  (22).  The  parameter 
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/fe 

\ dx 


)b 


The 


ra  is  evaluated  locally  by  using  the  definition  that  m 

resulting  values  of  a and  b are  used  in  Eqs.  (18)  and  (19)  to  detarmine 
6 and  A.  Ihese  results  are  then  used  In  Eqs.  (10)-(15)  to  calculate 
approximate  values  of  the  boundary-layer  parameters. 


STABILITY  ANALYSIS 

There  are  two  possibilities  of  approximately  relating  the  criti- 

v.cl  Reynolds  number  to  a boundary- layer  parameter.  The  first  is  the 

(3) 

Schlichting-Ulricb  correlation. 


= 645  exp (0.6a)  . 


(?3) 


(4) 

The  second  is  the  Lin  approximation. 


a = u’  c 

w 

R c^  = 25u'  , 

c w 


(24) 


(25) 

(26) 


where  the  velocity  has  been  normalized  with  the  subscript  prime 
denotes  d/(d  ^) , and  the  Reynolds  number  1 based  on  displacement 
thickness,  free-streara  velocity,  and  free-s*^ re.i»  viscosity. 

The  Schlichting-Ulrich  approximation  is  the  simplest,  and  it  will 
be  the  first  to  be  considered.  However,  a relationship  between  a and 
H must  be  determined  because  the  parameter  most  used  in  stability 
correlations  is  the  shape  factor  H.  From  Sqs.  (13)  and  (14)  one  can 
show  that  the  parameter  H is  related  to  a in  the  following  way: 


2a 


3_ 

5H 


(4H  + 21) 


- (4H  t 21)] 


288 

5H 


(2H  - 7) 


(27) 
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Note  that  the  critical  Reynolds  number  given  by  Eq,  (23)  depends  only 
on  the  shai'e  factor  H.  Therefore , n univeiL’a!  curve  for  as  a func 
tion  of  H can  be  calculated  that  includes  the  effects  of  pressure 
gradient,  variable  viscosity,  and  suction. 

]^et  us  investigate  the  Lin  approximation  to  see  if  the  observa- 
tion holds  true  there.  Using  the  profile  giveo  by  Eq.  (8),  one  can 
show  that  Eqs.  (2A)-(26)  reduce  to 


-T7a 


0.58 


(28) 


a = ac  , 


(29) 


R 


25a 
4 ’ 


(30) 


where  f is  defined  by  Eq.  (8).  Equation  (28)  could  be  used  to  deter- 
mine the  disturbance  wave  speed  c.  This  equation  states  that  the 
critical  layer  is  determined  only  by  the  profile  parameter  a,  which 
implies  that  c is  only  a function  of  the  shape  factor  H.  Equation 
(29)  implies  that  the  critical  Ct  wave  number  is  a function  of  H,  and 
Eq.  (30)  implies  that  the  critical  Reynolds  number  is  also  a function 

of  H alone. 

Two  different  approximations  of  the  critical  perimeters  for 
boundary-layer  stability  have  been  used  to  show  that  the  critical 
parameters  are  a function  oi  only  the  shape  factor  H.  This  conclu- 
sion is  independent  of  the  local  similarity  assumptlcn  of  the  previous 
section.  The  implication  for  correlating  experimental  ana  iasign  data 
is  encouraging.  For  example.  Fig.  1 (taken  from  Ref.  4)  sh-ws  a cor- 
relation of  R for  various  pressure  gradients  and  suction  effects  with 
the  shape  factor  H.  This  correlation  and  Eq.  (23)  are  in  fair  agree- 
ment when  2.3  ^ H s 2.70. 

There  are  indications  that  h s not  the  only  possible  universal^ 
parameter.  A parameter  discussed  and  defined  as  u"(0)  - ^ 
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in  Ref.  2 has  been  shown  to  be  an  effective  correlation  parameter.  In 
the  present  analysis,  the  relationship  between  u"(0)  and  a is 

u"(0)  = 6(2  - a)  . (31) 

This  indicates  that  u"(0)  is  related  to  H through  Eq.  (27).  In  compar- 
ing Eq.  (27)  with  Eq.  (31),  one  may  conclude  that  u"(0)  is  the  preferred 
parameter  for  correlating  theoretical  results  derived  from  integral 
solutions  because  Eq.  (31)  would  result  in  a simpler  computation.  How- 
ever, there  is  no  preference  in  correlating  theoretical  results  derived 
from  "exact"  numerical  computer  solutions  because  both  parameters  are 
equally  easy  to  calculate.  For  correlating  experimental  results,  one 
may  prefer  H as  the  parameter  because  it  is  easier  to  deduce  from 
measurements  than  u"(0). 

Another  interesting  consideration  that  can  be  seen  from  Eq.  (23) 
is  that  small  changes  in  the  profile  parameter,  a.  can  produce  consider- 
able changes  in  the  critical  Reynolds  number.  Changes  in  a are  ex- 
ponentially amplified  to  produce  larger  changes  in  R^.  Further,  it  can 
be  shown  that  changes  in  wall  heating  and  suction  have  a larger  effect 
than  small  changes  in  pressure  gradient.  It  is  not  apparent  that  the 
Lin  approximation  would  substantiate  this  conclusion  to  the  extent  of 
making  it  cogent.  However,  it  can  be  shown  that  the  Lin  approximatiop 
would  produce  simil-”-  results,  small  changes  in  heating  and  suction 
produce  larger  changes  in  R^,  and  the  numerical  results  of  Ref  2 have 
similar  indications. 

The  reason  for  these  conclusions  is  that  the  derivative  of  R^  with 
respect  to  either  the  wall-heating  parameter  or  the  suction  parameter 
is  larger  than  the  derivative  wich  respect  to  the  Pohlhausen  parameter 
by  the  factor  a,  ti:e  profile  parameter.  The  additional  implication  is 
that  the  larger  the  parameter  a or  the  more  favorable  the  pressure 
gradient,  the  more  stabilizing  are  the  effects  of  wall  heating  and  suc- 
tion. This  trend  is  in  agreement  with  one  established  with  numerical 
results  in  Ref.  5. 

It  is  interesting  to  compare  the  relative  effects  of  heating  and 
suction.  The  criterion  for  the  comparison  is  how  large  must  the  physical 
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boundary  conditions,  and  pv^/pU^,  be  so  that  both  ())  and  Q in  Eq. 

(12)  have  an  order  of  unity  magnitude.  According  to  Refs.  6 and  7, 
surface  conditions  of  this  magnitude  have  large  effects  on  the  boundary- 
layer  flow.  For  an  order  unity  <j>  and  = 60"F,  the  surface  has^to  be 
heated  to  at  least  lbO°F.  For  an  order  unity  Q and  pU^x/U^  ~ 10  , the 
suction  parameter  Pv^pU^  has  to  be  Pv^/pU^  ~ 10  Thus  suction 
appears  to  be  more  effective  than  heating  in  stabilizing  the  boundary 
layer.  The  Implication  here  is  that  it  is  much  more  difficult  to  pro- 
vide a system  that  will  produce  wall  temperatures  on  the  order  of  160“F 
than  to  provide  a suction  velocity  that  is  a tenth  of  a percent  of  the 
free-stream  velocity. 

DESIGN  IMPLICATIONS 

For  an  arbitrary  body  shape,  the  results  presented  here  can  be 
used  to  indicate  the  critical  boundary  layer  stability  parameters. 

For  example,  the  Pohlhausen  parameter  X can  be  calculated  using  the 
procedure  discussed  in  a previous  section.  Using  results  similar  to 
that  given  in  Fig.  2,  one  can  determine  the  profile  parameter  a for  all 
heating,  <|)  ?*  0,  and  arbitra’-y  suction  Q 0.  The  value  of  a is  used 
in  Fig.  3,  which  is  a plot  of  Eq.  (27),  to  determine  a value  of  H;  this 
H is  used  in  Fig.  1 to  determine  a value  of  R^.  Note  that  the  trends 
given  by  the  approximate  analysis  are  correct.  Increasing  wall  heat- 
ing, one  will  find  from  Fig.  2 that  a increases  for  a given  X.  Figure 
3 shows  that  as  a increases  H decreases,  and  Fig.  1 shows  that  as  H 

decreases  R increases, 
c 
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III.  CONCLUDING  REMARKS 


The  effects  of  wall  heating  and  ction  on  the  stability  of  a two- 
dimensional,  laminar  boundary  layer  have  been  studied  by  using  an  in- 
tegral analysis  of  the  boundary  layer,  and  the  Schlichting-Ulrich  and 
Lin  approximations  to  determine  rhe  critical  Reynolds  number.  It  was 
shown  that  the  effects  of  pressure  gradient,  wall  heating,  and  suction 
could  be  incorporated  into  a universal  curve  of  the  critical  Reynolds 
number  as  a function  of  either  shape  factor  H or  the  parameter  proposed 
by  Wazzan,  u"(0).  Further,  it  was  indicated  that  small  changes  in  wall 
heating  and  suction  produce  large  changes  in  the  critical  Reynolds  num- 
ber and  that  suction  is  more  effective  in  stabilizing  a boundary-layer 
flow  than  wall  heating.  It  was  also  shown  that  changes  in  the  critical 
Reynolds  number  produced  by  changes  in  pressure  gradient  were  smaller 
than  those  produced  by  changes  in  wall  heating  and  suction.  The  acaump- 
rion  of  local-similarity  was  used  to  discuss  a met.nod  for  locating  the 
critical  Reynolds  number. 

Although  the  parameter  crucial  to  designers  is  the  location  of 
transition  and  not  so  much  the  critical  Reynolds  number,  it  is  possible 
that  the  above  result  can  be  useful  in  correlating  transition.  A fur- 
tiier  study  of  existing  data  and  correlations  is  required. 
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